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In this work we investigate the effects of primordial magnetic fields on cosmic microwave back- 
ground anisotropies (CMB). Based on cosmological magneto-hydro dynamic (MHD) simulations Q 
we calculate the CMB anisotropy spectra and polarization induced by fluid fluctuations (Alfven 
modes) generated by primordial magnetic fields. The strongest effect on the CMB spectra comes 
from the transition epoch from a turbulent regime to a viscous regime. The balance between 
magnetic and kinetic energy until the onset of the viscous regime provides a one to one relation 
between the comoving coherence length L and the comoving magnetic field strength B, such as 
L ~ 30 (B/l 0~ 9 Gauss) 3 pc. The resulting CMB temperature and polarization anisotropies for the 
■ initial power law index of the magnetic fields n > 3/2 are somewhat different from the ones previously 

y—{ ' obtained by using linear perturbation theory. In particular, differences can appear on intermediate 

scales I < 2000 and small scales I > 20000. On scales I < 2000 the CMB anisotropy and polarization 
. spectra are flat in the case of our nonlinear calculations whereas the spectra have a blue index 

^ ' calculated with linear perturbation theory if we assume the velocity fields of baryons induced by the 

magnetic fields achieved Alfven velocity due to the turbulent motions on large scales in the early 
universe. Our calculation gives a constraint on the magnetic field strength in the intermediate scale 
CN ■ of CMB observations. Upper limits are set by WMAP and BOOMERANG results for comoving 

magnetic field strength of B < 28nGauss with a comoving coherence length of L > 0.7 Mpc for the 
most extreme case, or B < 30 nGauss and L > 0.8 Mpc for the most conservative case. We may 
also expect higher signals on large scales of the polarization spectra compared to linear calculations. 
The signal may even exceed the B-mode polarization from gravitational lensing depending on the 

pH . strength of the primordial magnetic fields. On very small scales, the diffusion damping scale of non- 

linear calculations turns out to be much smaller than the one of linear calculations if the comoving 
magnetic field strength B > 16 nGauss. If the magnetic field strength is smaller, the diffusion scales 
become smaller too. Therefore we expect to have both, temperature and polarization anisotropies, 

+a • even beyond I > 10000 regardless of the strength of the magnetic fields. The peak values of the 

temperature anisotropy and the B-mode polarization spectra are approximately 40 fiK and a few 

. . ' /xK, respectively. 



I. INTRODUCTION 



Magnetic fields have been observed in many galaxies and galaxy clusters |2|, |3| ■ Observations have revealed that 
these magnetic fields typically have a few /zGauss strengths and relatively large coherent scales, i.e., a few tens of 
kpc for clusters of galaxies and a few kpc for galaxies. It is one of the great challenges for modern astronomy to 
understand the origin of these magnetic fields. 

Perhaps the most conventional scenario of generating such magnetic fields is as follows. First, small seeds of the 
magnetic fields are produced due to the Biermann battery mechanism. Although the resultant magnetic fields are very 
weak, those are amplified by the dynamo process £| (for a comprehensive review see Eventually these magnetic 
fields are spread by Supernova winds or AGN jets into inter-galactic medium. 

However, the fact that magnetic fields with very large coherent scales are observed in galaxy clusters or high redshift 
galaxies with a strength of a few ^Gauss casts some questions on this standard dynamo scenario 0, El : how can these 
coherent magnetic fields be spread into inter-cluster medium, and is it possible for the dynamo amplification to take 
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place fast enough in high redshift galaxies? An alternative possibility to the dynamo scenario is the generation 
of magnetic fields in the early universe. There are in fact many previous works which suggest the generation of the 
magnetic fields in the early universe, e.g. during an inflation or at cosmological phase transitions (QCD or electroweak 
phase transition). For a detailed review, see [3. 

If magnetic fields are generated in the early universe, we need to understand their evolution in the expanding hot 
universe to obtain their field strength and structure at present. Several studies of evolution of magnetic fields have 
been performed by employing the linearized equations of magneto-hydro dynamics (MHD) These studies found 

that the growth of magnetic fields in the early universe is rather complicated such that the cosmological viscosity plays 
an important role M. Magnetic field energy is dissipated by the viscosity due to neutrinos and photons, equivalent 
to Silk damping [l(| for density fluctuations of baryon-photon fluid. However, the damping efficiency is different for 
different MHD modes. The Alfven and the slow modes are damped less than the fast modes so that the Alfven and 
the slow modes survive at small scales where the dissipation is effective. Meanwhile, it is also expected that the 
nonlinear effects contribute to the evolution of magnetic fields since equipartition between the magnetic field and the 
fluid will be established and the magnetic fields will cascade from large scales to small scales. Nonlinear effects on 
Alfven modes in the presence of viscosity is investigated analytically for particular configurations [£j and only a little 
damping of Alfven modes is found in these situations. 

Recently, Banerjee and Jedamzik 0, studied the evolution of Alfven modes using MHD simulation in the 
expanding universe including dissipation due to diffusion and neutrino/photon-drag due to free streaming neutral 
particles. The evolution of the primordial magnetic fields is solely determined by the kinetic Reynolds number, which 
is defined by R = v 2 /Lf where v is the fluid velocity, L is the length scale and / is the fluid dissipation (which 
is discussed in more detail in the following section), can be divided into three different regimes: turbulent regime 
(R 3> 1), viscous (diffusion) regime and (viscous) free streaming regime (R < 1). Furthermore, analytic expressions 
for the growth of the magnetic coherence length due to small wavelength damping and helicity conservation were 
found. Under optimistic assumptions it is possible for magnetic fields produced during a QCD phase transition with 
an initial coherence length of 1 pc to have attained a kpc coherence length at present. 

The existence of the primordial magnetic fields leaves traces on various cosmological phenomena. Investigating the 
cosmological effects of the magnetic fields in the early universe, i.e., the effects on Big Bang nucleosynthesis (BBN), 
structure formation and cosmological microwave background (CMB) anisotropies, we can constrain the strength and 
the structure of primordial magnetic fields. 

Among them, the traces of primordial magnetic fields on CMB anisotropies and polarization are of particular 
interest since they provide information of not only the magnetic field energy but also the coherence length. On the 
contrary, BBN gives only a limit on the magnetic field energy density (B < 7 x 10~ 5 Gauss at present) jl2l |. 

The coherence length is an especially important observational quantity if one wishes to understand the origin of 
cosmic magnetic fields because the length scale strongly correlates with the production mechanism. The primordial 
magnetic fields produced by a causal mechanism during, for example, a phase transition are limited by the length 
scale which corresponds to the horizon scale at the epoch of the magnetic filed generation. In contrast, the magnetic 
fields produced in the inflation epoch are expected to have the scale-invariant spectrum such that the amplitudes of 
the magnetic fields at the horizon crossing are the same over the all scales. 

There have been already several works to set constraints on the primordial magnetic fields by using current CMB 
data. But so far, non of these considered the full non-linear evolution of Alfven modes. Barrow et al. set 
a limit on the homogeneous (coherent) magnetic field by using COBE data. Such a homogeneous magnetic field 
produces large scale anisotropic pressures and these pressures require an anisotropic gravitational field to support 
them. This anisotropic gravitational field should generate anisotropies in CMB temperature. They obtain a limit of 
B < 10 _9 Gauss 0| at present. 

Primordial magnetic fields with a coherence length comparable to the horizon size at the last scattering surface 
(LSS) affect the sound speed of the baryon-photon fluid and change their acoustic oscillations. This effect might be 
observed as the modification of the acoustic peaks in the CMB angular power spectrum. Adams et al. found that 
this modification should be detectable by WMAP and PLANCK if the today's magnetic fields are B > 10 _8 Gauss 
with the coherence length larger than the horizon at recombination. 

A distortion of the CMB energy spectrum from the black body shape also gives a constraint on primordial magnetic 
fields. The magnetic fields with a small coherent scale are dissipated due to viscosity of the baryon-photon fluid and 
this dissipation energy distorts the black body spectrum of CMB [13 . The spectrum distortion is described by either a 
chemical potential /jora Compton y-parameter. If the dissipation process takes place in the very early universe when 
thermal equilibrium is maintained, it does not affect the CMB spectrum. As the temperature of the universe goes 
down, however, thermal equilibrium can no longer be maintained. Instead, kinetic equilibrium is preserved since only 
the Compton scattering process which conserves photon number is effective. If dissipation occurs during this stage, 
we expect to have fj, distortion. As the temperature drops more, eventually even kinetic equilibrium can no longer 
be maintained and we expect to have y distortion due to the dissipation. COBE-FIRAS gives the CMB chemical 
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potential constraint of /i < 9 x 10~ 5 . This limit corresponds to the magnetic field constraint of B < 3 x 10 _8 Gauss at 
present on a comoving scale of 400pc at redshift z>2x 10 6 On the other hand, the constraint by the Compton 
y-parameter, which is y < 1-5 X 10~ 5 , leads to the constraint on the magnetic fields of B < 3 x 10~ 8 Gauss on 600kpc. 

Another important constraint can be given from small scale CMB temperature anisotropics and polarization. Mag- 
netic fields induce peculiar velocities in photon-baryon fluids driven by Alfven modes. These peculiar velocities gener- 
ate temperature anisotropics and polarization in CMB. Observations by WMAP already constrain these anisotropies 
induced by Alfven modes to be smaller than the primary anisotropies due to acoustic oscillations of the photon- 
baryon fluid. Therefore, one might conclude that there is no chance to measure magnetic field sourced temperature 
fluctuations. However this is not the case for small scale anisotropies. Primary CMB anisotropies and polarization 
at the recombination epoch are expected to have a lack of small scale power due to the diffusion damping of the 
photon-baryon fluid (Silk damping). On the other hand, calculations of linear evolution of Alfven modes show less 
damping [la. ITtI ITsl fl{^ | . These Alfven modes generate CMB temperature anisotropies and polarization. Therefore 
we can likely expect to have dominant contribution on small scale CMB anisotropies from the magnetic fields if they 
exist. Using the WMAP result and taking into account both scalar and Alfven modes, one can constrain primordial 
magnetic fields to be B < 3.9 x 10 _9 Gauss on IMpc at present by using the Markov Chain Monte Carlo method [l^ . 




For the same reason, polarization induced by Alfven modes from primordial magnetic fields may dominate on small 
scales 20] . It is known that the polarization can be decomposed into two modes, i.e., the E-mode polarization (gradient 
component) and the B-mode polarization (rotation component) |2l|. Observationally they are distinguishable. It is 
very interesting that Alfven modes only produce B-mode polarization since its perturbations are of vector type. On 
the contrary, scalar type perturbations, which form the large scale structure of the universe, only produce E-mode 
polarization. It has been known that the gravitational lensing effect on the primary E-mode polarization caused by 
structure of the universe produces B-mode polarization on intermediate scales (a few tens of arc minute angular size) . 
However, polarization of the magnetic field origin can be dominant as a B-mode on small sales. 

A few authors have investigated the generation of B-mode polarization either by an analytic treatment with em- 
ploying the tight coupling approximation 22] or by direct numerical calculations |l3, 13 • F° r the evolution of Alfven 
mode perturbations, however, they all apply linear analysis. In this paper, we study the evolution of Alfven modes 
using recent MHD simulation by Banerjee and Jedamzik |ll 1 1 11] - which allows us to incorporate the fully non-linear 
and self consistent spectra of the Alfven modes to calculate resulting CMB temperature anisotropy and polarization 
spectra. 

This paper is organized as follows. In Sec. II, we summarize the evolution of Alfven modes based on the numerical 
MHD simulation of Banerjee and Jedamzik. We separate the evolution into three regimes, i.e., turbulent, viscous 
and free streaming regimes. In each regime, the evolution of velocity fields is carefully investigated. In Sec. Ill, 
the power spectrum of the velocity fields is calculated. Using this power spectrum, in Sec. IV, we compute CMB 
temperature anisotropy and polarization spectra. Sec. V is devoted to discussion, in which we compare our results to 
those obtained with linear perturbation theory. We give our conclusions in Sec. VI. Throughout the paper, we take 
WMAP values for the cosmological parameters, i.e., h = 0.71 (H = hx lOOKm/s • Mpc), T = 2.725K, fl h h 2 = 0.0224 



Our final goal is to investigate the effect of the primordial magnetic fields on the CMB temperature anisotropics 
and polarization. Magnetic stresses produce vortical modes, such as the "Alfven mode" and "slow magnetosonic 
mode" in the ionized fluids. As both are very similar we will refer to them as Alfven modes, henceforth. These Alfven 
modes generate additional temperature anisotropies and polarization by Doppler shift. Therefore, we need to know the 
evolution of Alfven modes first. Unlike previous works which employed a linear perturbation approximation [l6l . ll7l.l22] ] , 
Banerjee and Jedamzik [Hlll| recently investigated the nonlinear evolution of magnetic fields and Alfven modes using 
numerical MHD simulation. These simulations cover the three different damping regimes appearing in the early 
universe: turbulence, viscous, and free streaming. Based on their work, we summarize the evolution of the velocity 
fluctuations (Alfven modes) of the ionized fluid in this section. 

The growth of Alfven modes is affected by the interaction with photons and/or neutrinos in the early universe. 
Here, we consider the interaction with photons only as neutrinos are already decoupled from the cosmic evolution at 
the epoch of LSS. 

The MHD equations in the expanding universe including diffusion due to the photon background are given by [| 




and Q M h 2 = 0.135 [23| . 



II. EVOLUTION OF ALFVEN MODES 




(1) 
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(9B 1 

— + 2HB = -V x (v x B), (2) 
at a 

where pb is the baryon density, pt, Pt and c s t denote the total density, pressure and sound velocity of photons and 
baryons, H is the Hubble parameter and f is the dissipation term. The velocity v of photon-baryon fluids here is the 
Alfven mode, and as such is incompressible so that v satisfies V • v = 0. The dissipation term, f , is written as 

f = a- 2 KP + P)" 1 (V 2 v), (3) 



where the shear viscosity v is given by |24 
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V = i5^ Lmf P' ( 4 ) 



Here p 7 is the radiation energy density and L m f p is the photon mean free path, L m f p = (or^e) -1 where ctt is the 
Thomson cross section and n e is the free electron density. 
By using V • B = 0, Eq. is rewritten as 

_^ + _( v • V)v + (1 - 3c 2 st )Hv - -(v A ■ V)v A = f, (5) 
at a a 



where va is the Alfven velocity, va = B/ Ak /47r(p t + 

The dissipation term, f which assumes multiple scattering between photons and baryons is valid only on scales 
larger than L m f p . Since L m f p evolves rapidly as a 3 , it soon exceeds the wave length of a particular mode we consider. 
In other word, the comoving wave number of the velocity, k, becomes a/k < £ m fp- Subsequently, we need to treat 
baryons (ionized fluids) and photons separately. 

Once the diffusive description becomes invalid, Alfven modes will be damped due to free streaming background 
photons. Note that the damping time on the integral scale in the viscous regime is longer than the Hubble time 
and dynamic evolution of the magnetic field is 'stalled'. In the free streaming regime, the dissipation process can be 
described as radiation drag whose force is proportional to the fluid velocity v. The coefficient a of the drag force is 
given by 25] 



mfp- (6) 



3 pb 

Therefore, the Euler equation for the baryon fluid with free streaming photons is given by 



^ + - (v • V)v + Hw - - (El±Pl ) (v A • V)v A = -av. ( 7 ) 



dt a a \ Ph 

Let us now consider a characteristic scale on which most of the magnetic energy exists. This characteristic scale 
corresponds to the peak of the magnetic field power spectrum and can be defined by using the two point correlation 
function £(r) = (|B(x + r)B(x)|) as 



1 

m Jo 



dr£(r). (8) 



We refer to this scale as the integral scale and define the comoving wave number of the integral scale as h ln t = a/Li nt - 
The integral scale grows in time. In the turbulent regime a direct cascade due to non-linear interactions transfers 
energy from the integral scale to the much smaller viscous scale, where it is lost to heat. In the free streaming regime 
the integral scale may also grow due to the dissipation of flows on the integral scale itself. In both regimes the damping 
time scale is given by 

L a 

teddy ~ — ~ -T- ■ (9) 

V kv 

By comparing this eddy time scale on the integral scale with the cosmic time H~ l , we can judge whether non-linear 
cascade processes are important or not. If once tJJddy = AntA> < B r 1 is satisfied, magnetic energy dissipation becomes 
effective and magnetic energy on scales L < L mt is lost. 

Now let us classify the evolution of the MHD nonlinear Alfven modes into three regimes, i.e., turbulent, viscous, 
and free streaming regimes, following the results of the numerical simulation. In the following argument, we mostly 
consider the evolution of Alfven modes at the integral scale. 
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A. Turbulent regime 

In the early universe, since the photon mean free path L m f p , which is proportional to l/n e , is very short, the effect 
of viscosity on Alfven modes at the integral scale is negligible until the mean free path becomes large enough to have 
Reynolds number R ~ 1 on L lnt . For R > 1 we can ignore the dissipation term f in Eq. (J5J. The advective term 
(va • V)va/o drives the fluid velocity in this regime, which we call the turbulent regime. 

Ignoring the cosmological expansion term, we find that the fluid velocity eventually approaches an equipartition 
state: 

v = v A . (10) 

This behavior is consistent with the MHD numerical simulation P, [n| . Note that it is not clear whether the fluid 
velocity approaches an equipartition state on the scales larger than the integral scale or not since it takes longer on 
larger scales. We will discuss this point in Section 3. A. 

If magnetic fields do not decay, the Alfven velocity stays constant in time since B oc a -2 for the adiabatic expansion 
and pt oc a -4 during the radiation dominated epoch. Accordingly the eddy time scale at the integral scale evolves as 
'•eddy = ^int/^A oc Li n t/ \H / y/pt) oc a - On the other hand, the cosmic time 1/H is proportional to a 2 in the radiation 
dominated epoch. Therefore soon or later the cosmic time exceeds the eddy time, which allows for turbulent decay 
of the magnetic fields. 

Let us go into details about the evolution of the integral scale. By definition the integral scale corresponds to the 
peak location of the energy power spectrum of the magnetic field. Correspondingly, the Alfven velocity, which is 
proportional to the amplitude of the magnetic field, peaks at the integral scale. Since we assume a blue magnetic 
power spectrum, i.e n > 0, the eddy turnover time on the scales larger than the integral scale is longer than i^dy" 
Therefore the condition i e ddy = 1/H is first satisfied on the integral scale and then gradually moves towards larger 
scales. Meantime, the cascading decay of the short wavelength modes (k > k- m t) shifts the integral scale towards 
larger scales. 

The underlying physical process of the cascade decay is following. When nonlinear effects become prominent, the 
magnetic field energy and the fluid kinetic energy achieve equipartition. The flow eddies break into the smaller eddies. 
Hence the kinetic energy is transported from large scales to small scales by a nonlinear cascade (Kolmogorov process). 
The transported energy is changed into heat at the scale where the dissipation process is effective. Consequently, the 
magnetic field energy at the integral scale is converted into heat which is called direct cascade. 



B. Viscous regime 

As the universe evolves, the photon mean free path, L m f p oc a 3 , becomes larger. Dissipation due to photon 
interaction becomes efficient and velocity fluctuations are damped by the photon drag. Eventually the dissipation 
term dominates the advective term at the integral scale in Eq. ©. This is the second regime, which is refereed as 
the viscous regime. During the viscous regime, the eddy time at the integral scale is always larger than the cosmic 
time due to the decay of the fluid velocity v, which makes the eddy time longer. 

The transition epoch from turbulent to viscous regimes on the integral scale can be determined by comparing the 
advective term (v • V)v/a and the dissipation term f = (v/{pt + Pt)) V 2 v/a 2 . During the turbulent regime, the 
amplitude of the fluid velocity in the advective term is equal to the Alfven velocity va a s shown in the previous 
subsection. Therefore the advective term can be written by using the integral scale as (v • V)v/a ~ u A /i; n t. The 
dissipation term can be also rewritten as (y/(pt + p t )) va/ Lf nt . At the transition epoch, these two terms become equal, 
which yields 

U A = 7 — - — ry — = fH 1 - ( n ) 

[Pt + Pt)Mnt OL iat 

Here we assume the radiation domination, i.e., p t + Pt = 4p 7 /3. Until the beginning of the viscous regime, the eddy 
time teddy = L- mt /vA is equal to the cosmic time 1/H, which leads to va — \J H/ (n e aT)- Inserting the definition of 
lA into this equation, the transition redshift z t - v can be obtained as 

z t _ v ~ 6 x 10 6 BZ 9 2 ~ 6 x 10 6 ( ^ T ) , (12) 

9 V3.4 x 10 4 Mpc~V 

where _B_g is the comoving magnetic field normalized by 10 _9 Gauss, i.e., _B_g = (B/10 _9 )a 2 and 

fc int ~ 3.4 x lO^Z^Mpc^ 1 , (13) 
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is the comoving wave number of the integral scale at the transition epoch. Note that Eqs. l|12f) and l|13f) are only valid 
after e annihilation where £ m f p = l/a^n e , i.e., z < 10 8 . For earlier epochs Eq. I|ll[l must be evaluated numerically 
to find the redshift-B-field relations. 

In the viscous regime, we can ignore the advective term in Eq. (JSJ at the integral scale since the fluid velocity 
v decays due to the dissipation. We can also omit the expansion term since the evolution of the fluid velocity is 
controlled by the dissipation whose time scale is much faster than the cosmic expansion. Employing the terminal- 
velocity approximation, we obtain 



U A V U 



(14) 



Lint {Pt+Pt)L? lt 

Using the comoving wave number at the integral scale /c; n t , the fluid velocity can be described as 

2 

w=(pt+PtA^-. (15) 

V feint 

Hence it is found that the evolution of the fluid velocity becomes v cx a~ 2 . 

According to this solution, the eddy time at the integral scale is proportional to a 3 in both radiation and matter 
dominated epochs. Therefore the eddy time remains longer than the cosmic expansion time, 1/H which is proportional 
to a 2 , and a 3 / 2 in the radiation dominated epoch and the matter dominated epoch, respectively. As a result, no direct 
cascade occurs during the viscous regime and the integral scale does not grow. 



C. Free streaming regime 

The third regime is the free streaming regime, in which photon mean free path is larger than the integral scale. 
In this regime, photons and baryon fluids are decoupled and magnetic fields can amplify the fluid velocity. This 
amplification of the velocity makes the eddy time shorter until the eddy time becomes equal to the cosmic time 1/H. 
At this point the kinetic energy on L ln t is directly dissipated into heat and the integral scale shifts to larger scales. 

Let us first estimate the transition epoch from the viscous regime to the free streaming regime. The transition 
happens when the mean free path L m f p becomes equal to the integral scale Lint- Employing the comoving wave 
number of the integral scale in the viscous regime Eq. (|f 31) , we obtain the transition redshift as 

z v ^ { ~ 1.7 x 10 5 BZg /2 =s 1.7 x 10 5 f r ) , (16) 

9 V3.4 x 10 4 Mp C -V 

Using Eq. 1)15(1. and noting that a/feint ~ Lmfp during this epoch, we find that v has decayed to v ~ v\ at the 
transition. In the free streaming regime, we need to solve Eq. Q instead of Eq. JSJ for the evolution of the fluid 
velocity. Ignoring the cosmological expansion term, we obtain the solution with the terminal velocity approximation 

as 

v = jpt +Pt)h nt vi / 17 n 

Phaa 

The eddy time is longer than the cosmic time when the free streaming epoch begins since v is much smaller than 
va during the viscous regime due to dissipation. From the above solution, however, it is found that w cx a 2 in the free 
streaming epoch. Accordingly the eddy time evolves as a -1 . Therefore the eddy time soon becomes shorter than the 
cosmic time which is proportional to a 2 / 3 in the matter dominated universe. 

As is shown in the previous subsection, the integral scale does not evolve in the viscous regime. In the free streaming 
regime, the integral scale does also not change until the eddy time becomes equal to the cosmic time. 

When the eddy time is equal to the cosmic time on the integral scale, i.e. L- ln t/v — 1/H , the comoving wave number 
must satisfy the relation (using Eqs. (|17fl and iJBJ): 



K\ - 1, (is) 



where kg 1 is the comoving Silk scale defined as fcg 1 ~ y/L m f p /H/a. Note that in this regime, v^/ks is the comoving 
damping scale for Alfven modes. Unlike in the non-magnetized photon-baryon fluid overdamped Alfven modes survive 
below the Silk damping scale for weak magnetic fields as was first pointed out by @ and . While it is not clear either 
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the direct cascade or the diffusion process takes place first, the baryon velocity exponentially damps away bellow this 
scale. 

Substituting Eq. , which is the comoving integral wave number in the viscous regime and in the free streaming 
regime until the damping scale for Alfven modes, wa/^s, becomes equal to the integral scale, into Eq. I|18|) . we obtain 
the redshift below which further growth of the integral scale in the (viscous) free streaming regime occurs 

z = 1100 ( 7-7-77^ — ] =H00 — ^t) • (19) 
Vl6nGaussy VSMpc" 1 / 

Therefore if the integral scale once became larger than lOOkpc by the end of the turbulent regime, the integral scale 
did not change through viscous and free streaming regimes all the way until the recombination epoch. 

On the contrary, if the integral scale was smaller than lOOkpc, there was further growth of Li nt during the free 
streaming regime before recombination. The integral scale shifts to a larger scale in this case. The resultant integral 
scale at the recombination epoch depends on the slope of the magnetic field power spectrum. 

In Fig. 1 we show the areas of turbulent, viscous and free streaming regimes in the z — h m t plane. Here we employ 
Eqs. (|12fl . (|16fl . and l|19|) . Knowing the time evolution of the integral scale fc; n t, which depends on the magnetic 
field spectrum, we can draw an evolutionary track in this plane and easily understand when the transition between 
different regimes occurs. We will discuss this evolution in the next section. 




I 100000 10000 1000 

o 

redshift Z 

FIG. 1: The transition comoving scales as the function of the redshift. The dotted line is the transition line from the turbulent 
to the viscous regime, Eq. 1121 . The solid line is that from the viscous to the free streaming regime, Eq. 1161 . The dashed 
line represents the start of the direct cascade in the free streaming regime, Eq. 11911 . which is more or less when the diffusion 
process takes place. In the left side of this line there is no direct cascade or diffusion damping. In the right side the direct 
cascade or diffusion damping is occurring. 



III. ALFVEN MODE SPECTRUM 

In the previous section, we discussed the evolution history of magnetic fields and fluid velocity or Alfven modes. 
We separated the evolution into three regimes, i.e., turbulent, viscous, and free streaming. We analytically solved 
the MHD-Eulcr equation and acquired the relation between the fluid velocity v and the Alfven velocity va in each 
regime. Once v was obtained, we could estimate the eddy time t e ddy and compare it with the cosmic time in order to 
know whether field evolution occurred. 

However we can not describe the time evolution of the integral scale of the magnetic fields without knowing the 
magnetic field power spectrum. 

In this section, we find the power spectra of magnetic fields and the fluid velocities for a given initial shape of the 
magnetic field power spectrum. Being motivated by a causal mechanism of primordial magnetic field generation such 
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as QCD phase transition, we employ the power law spectrum with exponential cutoff for initial comoving magnetic 
fields as 

V%™(t u k) = ^(\B comov (t u k)\ 2 ) = B 2 comov (^fexp (- H , (20) 

where n is the initial power law index and fc c is the cutoff scale, and B comov is the amplitude of the initial comoving 
magnetic fields at the cutoff scale. Obviously, the cutoff scale corresponds to the integral scale at the initial epoch for 
a blue, n > spectra. 

To see the evolution of the fluid velocity induced by the magnetic field, it is convenient to define the power spectrum 
of the Alfven velocity Pa (A;) as 

v A (t,k) = ^(\ VA (t,k)\ 2 ) = ~ 1 4 n omov M)- (2i) 

2tt z 4-7T (p t + p t ) a 4 



A. Turbulent regime 



In the turbulent regime, we have shown that the magnetic field energy reaches an equipartition state with the fluid 
kinetic energy, i.e., v — va (Eq. l(TU|) '). 

Therefore we expect to have the following relation between the power spectrum of the fluid velocity V v (t, k) and 
the energy density of the magnetic field, i.e. the power spectrum of the Alfven velocity, V^t, fc), as 

P A (t,k) = P v (t,k). (22) 

Here the power spectrum of the fluid velocity is defined as 

-P v (k) = ^(\v\ 2 ). (23) 

Eq. (1221 is consistent with the results of the MHD simulation as is shown in FIG. [3 In this figure s is the time 
stamp and has the following relation with the redshift z, 

Si-- Rr> < M > 

where the subscript i denotes the initial value. Initially, the power spectrum of the fluid velocity is very different 
from the one of the magnetic field. However, the fluid velocity power spectrum soon catches up with the magnetic 
field one. Subsequently, they evolve in a very similar manner. In this figure, it is shown fluid velocity on large scales 
also immediately reaches the equipartition as well as the velocity on scales smaller than the integral scale. Therefore 
hereafter we assume the equipartition on all scales, i.e., Eq. (|22|) . Although it is immediate in terms of s, however, 
it may take longer time in physical time. If this is the case and the fluid velocity on larger scales has not caught up 
with the Alfven velocity by the end of the turbulent regime, we expect that the power spectrum of the fluid velocity 
V^k) has less power on large scales than the one of the magnetic fields V A {k) which is proportional to k n on large 
scales. The most extreme case is that turbulence does not work at all on large scales, in which we need to employ 
linear perturbation. Subramanian and Barrow [l6| obtain V v (k) oc k 5 for the linear calculation. Therefore the power 
law slope on large scales can be between n and 5. 

On scales larger than the integral scale, the power spectrum of the Alfven velocity keeps its initial amplitude. The 
power spectrum of the fluid velocity also soon reaches to this amplitude. 

Let us now estimate the evolution of the integral scale. In Sec. Ill Al we explained the reason why this scale shifts to 
larger scale with time during the turbulent regime. Now we would like to estimate the rate of this shift as a function 
of the power law index n. 

As we have shown in III Al the integral scale is the scale where i e ddy = 1/H is satisfied. In the turbulent regime 
Wdy = a/(kintv) = a>/{k- m tV a)- Therefore feint = aH/vA- The time evolution of the Alfven velocity at the integral 
scale, in the absence of dissipation, is estimated as va = B/sffh + Pt oc y/V A {kmt) x a°. On the scale larger than the 
previous integral scale, "Ps(fc) oc k n . Therefore the growth of the integral scale, when dissipation is also considered, 
is obtained as 



(25) 




1 2 5 10 20 50 

normalized comoving wave number 



FIG. 2: The evolution of the comoving energy power spectrum of the magnetic fields and the kinetic energy power spectrum 
of the fluids in the turbulent regime from the numerical simulation by 111 The x-axis is a comoving wave-number. The 
normalizations of x and y-axises are arbitrary. The thin and thick lines are the power spectra of magnetic fields and fluids, 
respectively. The dotted lines are the initial power spectra (r = 0.00) with the spectrum index n = 4. This figure shows that 
the energy densities become the equipartition, v — va after a few eddy time scale even on large scales. 



This analytic estimate is supported by the numerical simulation as is seen in FIG. [3 This figure also shows the direct 
cascade of the magnetic fields soon smears out the initial sharp damping in the power spectrum. 

The resultant power spectrum of the magnetic fields during the turbulent regime has power law shape on both large 
and small scales but with different indices. Matching the amplitude with the initial value on large scales, we obtain 

V A (t, k) = V A (t u fe) oc fe™, k<k int , (26) 
V A (t, k) = V A (U, fe ln t) [-;— I <x k m , k > fe int , (27) 



where m as observed in the numerical simulation is approximately m ~ —2/3 (for details see Q). 

Because of the equipartition, the evolution of the power spectrum of the fluid velocity is identical to the Alfven 
velocity as 



V v (t,k) = V A (t u k)ock n , fe< feint, (28) 
' _k 



V v (t, fe) = V A (t U fc in t) 1-^-1 CX fc"\ fe > feint. (29) 



B. Viscous regime 

When the relation v A = (p t + pt)vk/a is satisfied at the integral scale, the magnetic fields enter into the viscous 
regime. In the viscous regime the magnetic fields do not decay. Therefore the spectrum of the magnetic fields is frozen 
over the viscous regime and keeps its feature at the transition epoch from the turbulent to the viscous regimes t t - v as 

V A (t,k)=V A (t t - v ,k). (30) 

In the viscous regime no further growth of L- m t happens and the amplitude of the fluid velocity is proportional to 
square of the Alfven velocity amplitude as is shown in Eq. Ijl5|l . Note that Eq. (|15J) can be applied not only to the 
integral scale but also to all scales during the viscous regime. Therefore the power spectrum of the fluid velocity can 
be described using the power spectrum of the Alfven velocity. 

If the power law index n is smaller than 3/2, the contribution from larger fe (smaller scales) is negligible and we 
obtain 

t> n i\ ~ f((>t +Pt)a , A 2 / fc 2 ^- 1 ), fc<fcmt, r o-n 
7V*,fe)^ vk P A ftt-v,fe)j « | Jfc2(m _ 1)j k>k .^ (31) 
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If n is larger than 3/2, on the other hand, the contribution from the integral scale dominates in the convolution 
on scales larger than the integral scale. Then we can simply describe the power spectrum of the fluid velocity as (cf. 
Appendix A, Eq. I|A23|) and text below) 

V v (t,k) w ( ipt + Pt)a V A (t t - v ,k int )) (^-)ock, k< k int , (32) 

V fMnt / \ feint/ 



V v {t,k) « ^ Pt+ J k t)a v A {t t ^k)^\k^ 



k > feint. (33) 



Note that in Eq. (|32() . the slope may be steeper on large scales if the equipartition was not achieved during the 
turbulent regime. 

On scales smaller than the photon mean free path £ m f p , the diffusion approximation is no longer valid as is the 
case in the free streaming regime regardless of the value of the initial power law index n. Therefore we shall use the 
argument of Sec. Ill CI for the evolution of the fluid velocity. Employing Eq. JTJJ, we obtain 

V v (t,k) « ( {Pt+Pt)k V A (t t _ v ,k)) 2 oc fc^+D, fc > kt, (34) 

where fef = y/Ea/L m { p . For the definition of fcf, we put the factor V5 in order to match the power spectrum at kf. 
The estimated value of fcf at the recombination epoch is fcf = 1.2Mpc _1 (z/1100) 2 . 

C. Free streaming regime 

When the integral scale becomes smaller than the photon mean free path, i.e., fcf = fci n t, the free streaming regime 
begins. At the beginning of the free streaming regime, the eddy time at the integral scale is longer than the cosmic 
time. Therefore the integral scale does not change. If the integral scale is larger than lOOkpc, the eddy time is 
always longer than the cosmic time until recombination as is described in Eq. I|19|l . and further dissipation never 
occurs. In this case, the Alfven velocity v\ does not evolve and the fluid velocity is expressed in terms of v A as 
V = (pt + pt)kv\/ (phOta) (Eq. JTJJ). Accordingly the power spectrum of the fluid velocity can be described by the 
convolution of the power spectrum of the Alfven velocity. 

Following the argument we made in the previous subsection, we obtain the power spectrum of the fluid velocity as 

T> U U\~ ( (Pt+Pt)k , A 2 / k 2{n+1 \ fcf < fc < fcint, / qcn 

V ^ k ^{ phaa ?A(tv-f,*)j «| fc2(m+1)j k>k ^ (35) 

for the power law index n being smaller than 3/2. 
If n > 3/2, on the contrary, we obtain 

Vv{ t,k) « ( {Pt+Pt)hnt r A (tt- v ,k int )) (^-) ^5 kl<k<k . aU (36) 
V Phaa J V fcint/ 



V v (t,k) « ( {Pt+Pt)k V A (t v ^,k)Y 0Cfc 2 ^ +1 ), fc> fcint- CS7) 

On the scales larger than the free streaming scale, i.e., fc < fcf, the evolution of the fluid velocity still follows the 
solution of the viscous regime. Therefore 

V v (t,k) Pa ( {pt ^ Pt)a V A (t t . v , fcint)) (^-)o,k, fc<fc f , n > 3/2, (38) 



V v (t,k) » ~f^ pA ^- v ' fc )) <x k2(n ~ 1] > k<k U n<3/2. (39) 

Note that in Eqs. and l|S9l . the slopes may be steeper on large scales as is pointed out in the previous subsection. 

Finally, in the case of the integral scale at the transition epoch from the viscous to free streaming regimes being 
smaller than lOOkpc, the further processing of Li n t begins at the redshift of Eq. (|19|l . Once further decay takes place, 
the integral scale shifts to larger scales. Let us now estimate the evolution of the integral scale as a function of the 
initial power law index n. The integral scale can be written as fci„t = aH/v. The time evolution of the fluid velocity 
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at the integral scale is estimated as v = (pt + Pt)kv A / (pbaa) oc ki n tV A (ki nt )a 2 . On the scale larger than the previous 
integral scale, V A (k) oc k n . Therefore the time evolution of the integral scale is given by 

feint oc a~ 3 /(" +2 - ) in the radiation dominant, 

fcint oc a~ 5 /( 2 (™ +2 )) in the matter dominant. (40) 

We should note that H oc a~ 2 in the radiation dominated universe and H oc a~ 3 / 2 in the matter dominated universe 
while p t oc a~ 4 until the energy density of baryons dominate the one of photons, which happens much later than 
recombination. Therefore v A = B/^/pt + pt is still constant through this regime. This time evolution is consistent 
with the numerical result. 

Now we can estimate the power spectrum on scales where the direct cascade or diffusion damping during the free 
streaming regime occurs. Even if the cascade process occurs first, the power law slope of the decay immediately 
washes away by diffusion since the scales of the direct cascade and the diffusion are very close. Therefore we only 
consider the diffusion process here. The resultant power spectrum leads to 

Vv{t ,k) « ( ipt+Pt)hnt r A (k int ))\-^ k ^ 2 , k > k int . (4i) 

In the above power spectrum, the shape on the scales larger than the integral scale is identical to the one without 
further growth of the integral scale during free-streaming (cf. Eqs. (|35[1 - (jSHJ)- 



D. Summary of Velocity Power Spectra 

Let us summarize the fluid velocity power spectrum of Alfven modes at LSS for estimating power spectra of CMB 
anisotropics and polarization. There are three different cases corresponding to the integral scale on LSS. We refer 
them as "Case A", "Case B" and "Case C" for fe lnt < k { = 1.2Mpc~\ k f = 1.2Mpc _1 < k int < 8Mpc~\ and 
8Mpc _1 < fcint, respectively (cf. to FIG. H at z sa 1000). Here 1.2Mpc _1 and 8Mpc _1 correspond to amplitudes 
of comoving magnetic field strength as 30nGauss and 16nGauss. Case A is the case in which the viscous regime 
continues until recombination. For Case B, while transition from the viscous to the free streaming regimes occurs 
before recombination, the direct cascade process never happens on the integral scale. The diffusion process on the 
integral scale only takes place for Case C. 

Let us first list the power spectra with n > 3/2. 
Case A (Eqs. (23, E3> and <|31)l) 

V v (t,k) w f (Pt t Pt) V A (tt-v, feint)) flMocfc, fc<fe int , (42) 

V ^int / V fc int/ 

V v (t,k) w ( (Pt ^f )a 7M*t-v,fc)) txfc 2 ^- 1 ), k iat <k<k { , (43) 
V v {t,k) « ( {Pt+Pt)k V A {t t ^,k)]\k^\ k f <k. (44) 
Case B (Eqs. ® , © and 

V v (t,k) w f (Pt t Pt)< V A (tt-v,fent) > ) (^-)*k, k<k u (45) 



^rCint / \ &int , 

V v (t,k) * ( {Pt+Pt)kint V A (t t _ v ,k int )) (^-) ocfc 5 , fc f < fc < fc iat , (46) 
V Pbaa J V fci n t / 

V v {t,k) « ( iPt+Pt)k V A (t v _ u k))\k^\ k int <k. (47) 
Case C (Eqs. ® , and gTJ) 

V v (t,k) w ( iPt+ 1 Pt)a V A (t t . Vl k int )) (^-)ock, k<k u (48) 
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Vv (t,k) « f (Pt+Pt)fcint P A (tt-v,fcint)) (tM ocfe 5 , fc f <fe<fe int , (49) 
V Pb^a / \ hat J 

V v (t,k) * f (Pt+Pt)fcint P A (fc int )Ve- 2 ( fc /^) 2 , fc int <fc. (50) 

FIG. |21 shows the evolutionary track of each case. The Case A is represented as the dotted line. The Case B and 
Case C are the solid and the dashed lines. The slope of the tracks follow Eqs. I|25[) and 1)41) [I. Here we employ the 
spectral index of the initial magnetic fields n = 4. 

Next we list the power spectra with n < 3/2. 
Case A' (Eqs. (jSJ and ifSljl) 



7>„(t,fc) 
Case B' (Eqs. (GU) and 



pk P A (t t _ Vi fc)j «| fca(m -i) j 

2 



A; < feint, 

feint < fe < fef, 



V Pbcra / 



(51) 
(52) 



r v (t,k) 

V v {t,k) 



(^|^PA(*t-v,fe)) 2 OCfe 2 ^, 

\ 2 r fc 2(«+l) 
V A (ty-f,k) I OC I fc2 (m+l) 



vk 

(pt+Pt)k 



k < fef, 



Case C (Eqs. © (fef < fe < feint) and gTJ) 



fef < fe < feint, 
feint < fe- 



(53) 
(54) 



7Vt,fe) « )g ^A(*t-v,fc)^ ocfc 2 ^- 1 ), fc<fe f , (55) 

P v (t,k) « f (Pt+Pt)fc P A (t v -f,fe)) Ocfe 2 (" +1 ), fef <fc< fcmt, (56) 

^(t,fc) * f (Pt+Pt)fcint P A (fcint)Ve- 2 ( fc /^) 2 ; fcint<fc . (57) 

Note that on large scales (fe < feint), the power law slopes of above equations may be steeper as is pointed out 
before. In the most extreme case, P v (t, k) cx fe 5 instead of fc for n > 3/2, and V v (t, fe) oc k 2(n+ ^ instead of fe 2 '"- 1 ) for 
n < 3/2 as is expected from the linear analysis. 



IV. CALCULATION OF CMB TEMPERATURE ANISOTROPY AND POLARIZATION SPECTRA 



Alfven modes of baryon-electron fluid produced by the primordial magnetic fields generate CMB temperature 
fluctuations by a Doppler shift 0|. Moreover, the quadrupole component of the generated temperature anisotropies 
produces polarization due to Thomson scattering. It is known that CMB polarization can be decomposed into two 
parity independent modes, i.e., the E-mode (electric type) and the B-mode (magnetic type) |2l|. Among them, the 
B-mode polarization is not primarily produced by the scalar type perturbations which eventually form the structure of 
the universe and provide dominant contribution on the temperature anisotropies. Therefore the B-mode polarization 
possibly can be a good probe of primordial magnetic fields. 

First, we summarize the derivations of the temperature anisotropy and polarization spectra induced by Alfven modes 
with use of the linear perturbation theory \wi Il7l l22j . The Alfven mode in the linear MHD theory is identical to the 
vector mode of the cosmological fluid velocity in the cosmological perturbations. The vector metric perturbations in 
the linear perturbation theory can be written as (cf. (2(| also for notation) 



<5.9oi = -a 2 V t = -a 2 VQi, Sg^ = 2a 2 H Tij = 2a 2 H T Qij, 



(58) 
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FIG. 3: Evolution of the comoving integral scale of the magnetic field configurations with n = 4. The dotted line corresponds 
to Case A which undergoes the turbulent regime and the viscous regime before recombination. The solid line corresponds to 
Case B which passes through turbulent, viscous and free streaming regimes, with nevertheless no growth of the integral scale 
immediately before recombination occurring. The dashed line corresponds to Case C in which there is growth of L lnt in the free 
streaming regime shortly before recombination. The thin dotted line represents the transition line from the turbulent regime 
to the viscous regime. The thin solid line shows the evolution of kf, whereas the thin dashed line represents the start of growth 
of the integral scale in the free streaming regime. 

where V and Hx are the amplitudes of the vector mertric perturbations and Qi and Qij are the vector type mode 
functions which are defined as 



V 2 Q. t = -fc 2 Q 4 , V l Q, = 0, 



1 

2fc 



(ViQj + V 



3 3^1) 



(59) 



The amplitude «b of the vector components of the baryon fluid velocity v\,i can be represented with the mode functions 
as 



Vbi = v h Qi. 



(60) 



Employing the total angular momentum method [2]j, we obtain the CMB temperature anisotropics AT TT (l) induced 
by Alfven modes as 



AT 1 (I) = T 



l{l + l)C TT (l) 
2tt 



(61) 



2tt 



(62) 



(63) 



(64) 



21 + 1 



dij e 



r(tfc-VU (11) (*(»?-»/)) 



+ (fP (1) (k, rf) + ±kV\ (k( v - r/)) 



(65) 
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^Zt + f = [ dffe^tP^ (k, rj')^ (k( V - rf)), 



(66) 



E^\k,r)) = -5V6 / dr}'e- r WfpW(k,'q')£ l i ) {k{<r 1 -ri 1 )), 



(67) 
(68) 



ii (U) («) 



l(l + l) jl{ X ) (21) 
O „ 1 Jl W 



31(1 + 1) d fj t {x) 
2 dx 



(69) 



(x) 



32{x) ldh(x)_ 
x 2 x dx 



(70) 



dl 1) (x) = -^(l-l)(l + 2) 



3i( x ) 



(71) 



— exp [ — / rdrj 

v' 



(72) 



where ji{x) is the spherical Bessel function, 77 = J da/ a is the conformal time, subscript denotes the present epoch 
and f is the differential optical depth which is expressed as 

f = an e a T = a/L m { p . (73) 

For further notation the reader is referred to Ref. |27j . The conformal time r\ can be written in terms of redshift z as 

1 f dz 



77 = 



oJ yjl - tt M + fi M [l + (1 + z)(l + zc^l + 



(74) 



where z cq = 2.4 x 10 4 57m^ 2 is the redshift of the matter-radiation equality epoch. In the above expressions, we have 
taken two independent vector modes into account. The time evolution of i>b is given by the Euler equation. 

Combining these formulas with the linear perturbation theory, we can calculate temperature anisotropy and polar- 
ization spectra induced by the nonlinear Alfven modes by simply substituting the fluid velocity of Alfven modes v in 
Sec. [n]to v b — V of Eq. H65fl due to the Newtonian gauge condition |27J- Although, we obtain the Alfven modes from 
(fully non-linear) numerical simulation, we still can treat them as the vector modes because V • v = 0. 

For the recombination history, we assume the visibility function to be a Guassian as 



g(i]) = fexp 



fdrj ) = (27rc7 2 )- 1/2 exp 



(v - ?7lss) 



2 1 



(2a 2 ) 



(75) 



where 77LSS is the conformal time at the LSS and a is the width of the LSS. From the WMAP results, the redshift 
and the thickness of the LSS are 1 + z — 1089 and Az = 195 [2j|. This implies that a is 13Mpc for the cosmological 
parameters we use in this work. 

In our calculations the magnetic field power spectrum stays virtually unchanged after the transition from turbulence 
to viscous diffusion (i.e. for t > t t - v ), even in the case where k mt = lOMpc" 1 (Case C), since the free streaming regime 
before recombination is so short that there is hardly any time for Lj nt to grow. Since in this study we concentrate on 
causal spectra (n > 3/2) we mostly require the magnetic filed spectrum around the integral scale in order to derive 
the fluid velocity spectrum on all scales. We therefore take the numerical results in the fully developed turbulent 
regime at s — 2.49, which is expected to resemble that at tt- v closely. From these numerical spectra of the magnetic 
fields, we obtain the power spectra of Alfven modes using Eqs. I|44 |l -I|57 |l . We assume the equipartition between 
magnetic fields and fluid velocities during the turbulent regime, which makes the power spectra flat on large scales 
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FIG. 4: The temperature anisotropy power spectrum. The dashed line is the case of the magnetic fields with the comoving 
integral scale k lnt = l.OMpc" 1 with 32nGauss. The solid line is the case of feint = 1.5Mpc _1 with 28nGauss and the dotted 
line is the case of hint = lOMpc" 1 with 16nGauss. The comoving integral scale is related with the magnetic field strength by 
Eq. 1131 . The gray line is the linear result of the magnetic field strength 28nGauss at 1.5Mpc with spectral index n = 4. This 
line has exponential damping at smaller scales than the cut-off scale (for details see Sec. IV CV For purpose of comparison, we 
plot the temperature anisotropy power spectrum in the standard ACDM model (a thin line), computed using CAMB with the 
same cosmological parameters |2q|. 



(small Vs). The modification of the power spectra due to possible violation of this assumption will be discussed in 
the next section. 

For the calculation of polarization, we compute the source term PW by using the publicly available code CAMB 
[Lsl l28| in which we substitute Alfven modes obtained by the numerical simulation. 

Now we are ready to calculate CMB temperature anisotropy and polarization spectra. Hereafter, we fix the initial 
power law index n to be 4. Temperature anisotropy and B-mode polarization spectra are shown in FIGs. 01 and \5\ 
respectively. Here we calculate models with k[ n t = lOMpc -1 , 1.5Mpc -1 and l.OMpc -1 at the LSS, whose magnetic 
field strengths at the integral scale are 15nGauss, 28nGauss, and 32nGauss, respectively (see Eq. lO). Note that 
feint = lOMpc -1 corresponds to Case C, 1.5Mpc -1 to Case B and l.OMpc -1 to Case A in Sec. IIII Dl 

For purpose of comparison, we also plot the temperature anisotropies and polarization by using linear perturbations 
of the Alfven modes in FIGs. 0] and [SJ Here we adopt two types of the magnetic field strength. One is 32nGauss at 
l.OMpc -1 and another is 28nGauss at 1.5Mpc -1 . In both cases we take the spectral index as n = 4. 



V. DISCUSSION 



Let us investigate the features of the power spectra of CMB anisotropies and polarization obtained in the previous 
section. For the qualitative understanding of these spectra, we first develop analytic expressions of CMB anisotropies 
and polarization. 

Both spectra are induced from the power spectrum of the fluid velocity P v (k) at the LSS. In fact, employing the 
small angle approximation and considering the phase cancellation damping of anisotropies within the thickness of the 
LSS, we obtain [Hill 



(76) 

k=l/{rio— I7lss) 



AT BB {1) M T 



v^ 7 

2ka 



'p(i) 



(A:) 



(77) 



k=l/{r]o—iTLSs) 
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FIG. 5: The B-mode polarization power spectrum. The dashed line is the case of the magnetic fields with the comoving 
integral scale k lnt = l.OMpc -1 with 32nGauss The solid line is the case of k lnt = 1.5Mpc _1 with 28nGauss and the dotted line 
is the case of fent = lOMpc- 1 with 16nGauss. The comoving integral scale is related to the magnetic field strength by Eq. 
1131 . The gray line is the linear result of the magnetic field strength 28nGauss at 1.5Mpc with spectral index n — 4. Because 
this linear result uses the tight coupling approximation, it is plotted up to the ~ kj w 1.4 Mpc -1 . For purpose of comparison, 
we plot the polarization power spectrum in the standard ACDM model (a thin line), computed using CAMB with the same 
cosmological parameters. This anisotropies is produced by the gravitational lensing mainly. 

where V P (i)(k) is the power spectrum of the source term defined as 

P p(1) (fc)^(|pM| 2 ). (78) 

Here, l/fccr is due to the phase cancellation damping within the thickness of the LSS. 

In Sec. IHII we studied the velocity power spectra where we found that the shape and amplitude of the power spectra 
on the LSS depend on the regime of the comoving integral scale at recombination. Note that magnetic fields with the 
comoving integral scale fc; nt = l^Mpc" 1 at recombination are transiting from the viscous regime to the free streaming 
regime right at recombination as is shown in FIG. |3| At the recombination epoch, the integral scale is still in the 
viscous regime for k- mt < 1.2Mpc , while it has already been in the free streaming regime for h m t > 1.2Mpc . 

In the viscous or free streaming regimes, the velocity power spectra do not explicitly depend on the initial spectral 
index of the magnetic fields n if n > 3/2, which is the case we consider in this paper, as is shown in IIII Bl and IIII d 
(or summarized in IIII D|) . Even in this initially blue spectrum case, however, the initial spectral index still affects 
the growth rate of the integral scale (see Eqs. (|25l) for the turbulent regime and (|4U|I for the free streaming regime). 
But as is shown in FIG. |3| the integral scale does not change during the viscous regime and only changes little in 
the free streaming regime, provided the magnetic fields are sufficiently strong, i.e. the comoving B > 16 nGauss at 
recombination. Since the turbulent regime, where the integral scale mostly evolves, is long before the recombination 
epoch, we can ignore the growth of the integral scale and accordingly the dependence of the initial spectral index. 
Therefore we expect the velocity power spectra and the resultant temperature anisotropy and polarization spectra to 
show virtually no dependence on the initial spectral index n, for n > 3/2. 

A. Temperature anisotropy power spectrum 



Let us now discuss the temperature anisotropies. Following Sec. IIII Dl we investigate the temperature anisotropy 
spectrum for three separate cases corresponding to the integral scale on the LSS. 

First, we consider Case A in which the integral scale fc; n t < 1.2Mpc . The integral scale is still in the viscous 
regime at the recombination epoch in this case. Substituting the equations of Case A in Sec. 1111 Dl into Eq. I|76|). we 
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can derive the approximations of the temperature anisotropy spectrum produced by the magnetic fields as 

AT TT 



T 



104 



2ka V vk- mt I \k int 



fc=V( r /0-'?LSs) 



kit 



l.OMpc" 



-13/6 



\pK], l<hnt, 



(79) 



AT 



TT 



To 



104 



( (pt +Pt) a 



2ka 



l.OMpc" 



vk 

-13/6 



fe=i/('7o-»7LSs) 



-3/2 



\pK], l hlt <l<l f , 



(80) 



AT TT _ T V^_fiPt±Pt)_k 
V 2fccr \ pbaa 



V A (tt-v,k) 



31.5 



J, \ -1/6 / , \ m+1/2 



l.OMpc" 1 



fc = V(»)0-»?LSs) 

[fjK], 1 { <1, 



(81) 



where 



iint = feint (vo - Vlss) = 14500 



kit 



l.OMpc" 1 



(82) 



k = k{( Vo - Vlss) = 17800. (83) 

We find a very good agreement between these approximations and full numerical calculations as is shown in FIG. 
0] On scales I < l int and for ka > 1, Eq. l(79"|) suggests no I dependence of the temperature anisotropy spectrum, 
whereas for I < l int and ka < 1, although this region is not shown in FIG. 01 a residual small Z-dependence is expected 
to remain oc I 1 / 2 . Since m ~ —2/3, I dependence on scales I > If is also very weak. We expect to see some variation 
of the temperature anisotropy spectrum between k nt < I < If. The dashed line in FIG.^shows these features. 

Note that we assume the equipartition between magnetic fields and fluid velocities during the turbulent regime. On 
very large scales (small i's), we may expect violation of this assumption. Above the scale at which this assumption is 
no longer valid, the temperature anisotropy spectrum shows / dependence. In the most extreme case, the advective 
term can be completely ignored. Therefore the behavior of the power spectrum is described by the linear perturbation 
as is shown in Appendix (A29), and the temperature anisotropy spectrum is proportional to l^l 2 . In FIG. 4, this 
damping on small Vs (the blue spectrum) is shown for the linear case. For the nonlinear cases, we do not really know 
the location of the break from the flat spectrum and the power law index since they strongly depend on when and 
how the initial magnetic fields were formed. If there existed strong magnetic fields in the very early universe, the 
assumption of the equipartition is mostly valid, and the flat spectrum on small Z's is expected. 

The amplitude of temperature anisotropies in Case A is generally very large. Even in the case when k m t = 1.2Mpc _1 , 
which corresponds to the smallest amplitude of Case A (or transition between Case A and Case B), AT TT ~ 70/iK. 
This amplitude is much larger than the one measured by WMAP or BOOMERANG, i.e., AT TT ~ 40^K at I ~ 800 
[2^| or AT TT ~ 30/^K at I ~ 1100 [3(|, in which major parts of anisotropies should be explained as the primordial 
anisotropies. Only a small fraction of anisotropies can be interpreted as primordial magnetic field origin at most. 
Note that we ignore here the possible damping of the temperature spectrum on small Vs due to the non-equipartition 
between magnetic fields and fluid velocities during the turbulent regime. Accordingly the amplitude of temperature 
anisotropies at / = 800 and 1100 may be smaller than 70/iK. Even in the most extreme case with assuming the same 
power spectrum as the linear calculation, however, we still get AT TT ~ 20/iK at I = 1100, which is about 3/2 of the 
observed temperature anisotropies. Here we take the break scale l v is 1800 using the expression in Appendix, and 
employ the equation AT TT = 70(l/l v ) 5/2 fiK. 

Therefore we can safely rule out the possibility of Case A, i.e., ki nt < 1.2Mpc _1 and B > 30nGauss at the 
recombination epoch (see FIG.EJ) for n > 3/2. 
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FIG. 6: The magnetic field strength B-g for the comoving integral scale feint (Eq. I13H 'I. The gray region is the forbidden 
region by the WMAP results. Magnetic fields with comoving strength larger than 30nGauss at recombination are ruled out. 

Now we study Case B. In this case, the integral scale is 1.2Mpc _1 < feint ^ 8Mpc _1 and the temperature anisotropies 
on the LSS are given by the power spectrum in the free streaming regime. Substituting the equations of Case B into 
Eq. 0, we get 
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Again we find a good agreement between the above approximations and full numerical calculations, i.e., a solid line 
in FIG. 01 On the scales I < If and I > Zi nt , we expect no, or very little I dependence from Eqs. (|84() and (|86(l while 
the non-equipartition between magnetic fields and fluid velocities may change I dependence in the same manner as 
Case A. Between If < I < ^ n t, however, temperature fluctuations should increase as I 2 . In FIG.H the solid line stays 
flat on the scale I < 18000 ~ , and increases until I ~ 30000. However AT TT oc I 2 only on the scale I < 22000 ~ Z int 
and gradually decreases the gradient for larger /. 

The WMAP and BOOMERANG results have ruled out Case A. They also constrain Case B as fe int > 1.5Mpc~\ or 
B < 29nGauss at the recombination epoch if the equipartition assumption is valid. When the assumption is violated, 
we cannot set concrete limits for Case B. On the other hand, one might think that CMB anisotropies produced by the 
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mag netic fields can explain the small scale excess of the temperature power spectrum observed by the CBI experiment 
[31j. However, the power spectrum of CMB anisotropics we find is very flat on scales 2000 < I < 18000 even with 
considering possible violation of the equipartition. Therefore it is rather difficult to explain the CBI experiment which 
shows Silk damping at / = 2000 and increases the power to I = 2800. 

Finally let us investigate Case C, in which some further growth of the integral scale takes place before recombination. 
In this case, the power spectrum Eq. (|86(l is no longer valid due to the direct cascade while the rests of the power 
spectra Eq. I|84|) and l|85|) are the same. Therefore the difference between Case B and Case C appears in the power 
spectrum at scales smaller than the integral scale. Employing Eq. I|50|l . we obtain the temperature anisotropy 
spectrum as 
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TT 




ifk +Pt) fcint"" 
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The dotted line of FIG 01 corresponds to Case C. In this figure, we can find the temperature anisotropy spectrum 
stays constant for I < If and follows cx I 2 for I > If. The integral scale of this model is too large, i.e., /j n t = 1-5 x 10 5 , 
to see the difference with Case B. 

Note that the temperature anisotropy spectrum induced by the magnetic fields shown in FIG. ^ does not show any 
significant damping on large Vs (small scales) unlike the primary temperature anisotropy spectrum which suffers Silk 
damping at I ~ 2000. This is due to the fact that Alfven modes can survive even below the Silk damping scale •8]. 
The damping scale for Alfven modes is the integral scale. Our numerical simulations also show the damping below 
the integral scale is rather mild compared to the linear calculations. There exists the diffusion scale wa/^s below 
which the power spectrum is exponentially damped for Case C (and Cases A and B in extremely small scales), while 
we cannot see the damping within the range of FIG. 01 Further comparison with linear calculations will be made in 
Sec. WU 



B. B-mode polarization power spectrum 

We investigate the B-mode polarization spectrum shown in FIG. [21 For the quantitative understanding, we need to 
know the behavior of the source term PW defined in Eq. I|t)7[) . In our calculation, we employ CAMB to solve PW. 
However, it is rather easy to develop an analytic solution for the evolution of P^\ 

First of all, we divide the evolution into two stages. The first one is the period when the tight coupling approximation 
is valid. Once the wave length of the perturbations becomes shorter than the Silk scale, however, electrons and photons 
are decoupled and the tight coupling approximation is no longer valid. This gives rise to the second stage. 

In the first stage, the solution of the Boltzmann equations of polarization becomes = — (V6/4) in the 
limit of t — > oo (Eq. (89) of pTjp. The quadrupole (I — 2) component of the Boltzmann temperature hierarchy 

gives — (4V3/9) (k/r) Qj . Due to the tight coupling, electron (baryon) velocity follows the photon velocity as 

Oj 1 = v. Therefore the source term can be written as (Eq. (94) of [27j |^ 

P(D ^ l [e « - V6E^} = V } = — V 3 = ^-^Lv. (88) 
10 L 2 2 J 4 2 9 t 1 9a v; 

In the second stage, the E-mode component is damped due to diffusion. Therefore Q% induced by the Alfven 
modes (Eq. H65[) ~) v only contributes to the source term P^ in Eq. I|67|) . Namely, PW is induced by v. It is known 
that the oscillations of the temperature quadrupole, i.e., (k(v — v)) °f Eq. (|65Jl . suffers damping due to the phase 
cancellation within the optical depth. Note that the 1/Vka coefficient of Eqs. (176(1 and Q77JI appeared because of this 
effect. Therefore we expect to have a factor 1/ ^kL m f p /a in the source term P^\ Adding a numerical factor to fit 
with the simulation, we obtain 

PW « / . (89) 

lOA/fcLmfp/a 



which we refer as the decoupling approximation. 
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For the intuitive understanding, we plot the source term (k, ?7lss) of the model in which the magnetic fields 
have a scale-invariant power spectrum, i.e., k 3 P v (k) oc k°, with a magnetic field strength of nGauss in FIG.[7|by 
employing CAMB. In this figure, it is found that the numerical calculation follows the tight coupling approximation 
at scales larger than the Silk damping scale, fcs = ay/H/L m { p ~ 0.08Mpc _1 . It is also shown that the numerical 
calculation traces the decoupling approximation Eq. (|89() on scales smaller than the Silk damping scale. Gradual 
increase of the source term on scales k > kf = ^/5a/L m f p ~ 1.2Mpc _1 is due to the increase of the velocity v since 
the viscosity is no longer efficient and the fluid velocity can evolve on these scales (see Eq. I|17fll. 
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FIG. 7: The Thomson scattering source term Pi(k) at the recombination epoch obtained by CAMB for a fixed magnetic field 
with InGauss amplitude and the scale-invariant power spectrum (solid line) . The dashed and dotted lines are the tight coupling 
approximation, Eq. 1881 . and the decoupling approximation, Eq. 1891 . respectively. Note that the Silk damping scale here is 
fcs — 0.08Mpc -1 , which corresponds to the transition scale between tight coupling and decoupling and the free streaming scale 
is fc f ~ 1.2Mpc _1 . 

Now we are ready to calculate the B-mode polarization spectrum using Eq. (|77|l . since the source term P' 1 ' can be 
written by employing the fluid velocity v as Eqs. (|8*8l and Following our calculations of temperature anisotropies, 
we obtain a B-mode polarization spectrum for three separate cases, i.e., Cases A, B, and C whose velocity spectra 
can be seen in Sec. IIII Dl 

For Case A, the approximation of the B-mode polarization spectrum is expressed as 
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where 



For Case B, 



AT 



BB 



5.2 



"ant \ I L 



l.OMpc -1 



[fjK], l f < I, 



k = ksivo - Vhss) = 1800. 



Tn 



0.09 



Vtt / kL mtp \ f (p t +pt)a 



2ka V 3a 



-Va (*t-vj feint) i — 

^ feint / V feint 



feint \ ( I 



1.5MPC" 1 



\uK], I < h 

500 ' 



(93) 



(94) 



(95) 



AT 



BB 



26 



1 a ( (p t +pt) a 



2ka 100 fcimfp \ vk 



L.. \ -8/3 / , \ -1/2 

"-mt \ i 1 



"T^A^t-v, feint) 



,1.5Mpc~ 



k 

feint 



[pK], l s <l<h, 



fe='/( r /0-»?LSs) 



(96) 



AT 



SB 



ji / 3V^ 1 a ( (A +p t )fcint p ^ fc AV fc x r 
V 2fccr lOOfcimfp V p h aa t_v ' lnt / Vfe int 



9.9 



1.5Mpc~ 



-2/3 / ^ x 3/2 



— ) \flK], If < I < lint, 

''in 



(97) 



AT 



BB 



To 



9.9 



' 30? 1 a / (pt + Pt) fe 
2ka 100 feL m fp \ pbaa 



,1.5Mpc" 



-2/3 



-r A (t v -{,k) 



b m t < I- 



(98) 



Finally, for Case C, the difference with Case B only arises at the scales smaller than the integral scale as is the case 
of temperature anisotropics. Using Eq. (|50|l . we acquire the B-mode polarization spectrum as 
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Let us discuss our numerical results of B-mode power spectra in FIG. |3] The dashed, solid and dotted lines 
correspond to Cases A, B, and C. It is shown that all lines gradually increase until I ~ 4000, since AT BB oc I on 
scales / < Is for all cases. Note that we assume the equipartition between magnetic fields and fluid velocities during 
the turbulent regime. On very large scales (small Vs), we may expect violation of this assumption and steeper slopes 
for the B-mode power spectra below I < 2000 as is the case with temperature spectra. Case A shows continuous 
declining on I > 4000 which is consistent with our analytic estimate. The increase of AT BB on scales I > It = 17800 
for Cases B and C is caused by the efficient evolution of the velocity field in the free streaming regime due to the lack 
of dissipation term f. For both Cases B and C, AT BB decreases on scales I > l m t- We can see this decrease in FIG.|S] 
for Case B whose h n t is 22000. For Case C, we cannot find this trend because h n t = 1.4 x 10 s , which is beyond the 
range of FIG.|3 

From the temperature anisotropy spectrum, Case A and a part of Case B are excluded. However, the solid line in 
FIG03 which overshoots the B-mode produced by the gravitational lens effect, is not yet ruled out. Even the dotted 
line can provide dominant contribution as the B-mode polarization on scales I > 3000. 
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C. Comparison with Linear Calculations 

Here, we compare our results with linear calculations. Note that all previous works investigated the effects of the 
magnetic fields on CMB anisotropics and polarization are based on linear perturbation theory. 

In the Appendix, we summarize the linear perturbation analysis and the approximated estimations of temperature 
anisotropy and polarization spectra following the analysis by Subramanian and Barrow f|l6l l22|'). 

In FIG. 01 we plot the temperature power spectrum of the linear calculation with B = 28nGauss at k = 1.5Mpc _1 
and n = 4 as the gray line for comparison. It is shown that differences with the nonlinear calculation (the solid line) 
appear on scales I < 2000, and I > 18000. 

On large scales, we can ignore the diffusion in the early epoch. In the nonlinear calculation, v immediately 
approaches the Alfven velocity v = va due to the nonlinear coupling during the turbulent regime if there exist strong 
magnetic fields in the very early universe. On the other hand, the growth of the velocity induced by the magnetic 
pressure gradient in the linear calculation turns out to be much slower as is shown in the Appendix (Eq. 1|A9|0 . 
Accordingly we find a smaller temperature power spectrum (and polarization) for the linear calculation. Note that 
even for the nonlinear calculation, there may be the case in which the cquipartition could not be achieved by the end 
of the turbulent regime if the primordial magnetic fields on the large scale were not strong enough. In this case, we 
rather expect to have similar behaviors of the power spectrum to the linear one on large scales (small Vs). 

Once the viscosity becomes efficient, the velocity follows the linear solution even in the case of the nonlinear 
calculation (see Eqs. (|15|l and IjAlOJl ). Unlike the perturbations without the magnetic fields which suffer severe 
damping, the velocity (Alfven mode) induced by the magnetic fields can survive within the Silk scale for both nonlinear 
and linear calculations. The amplitude of the velocity is determined by the Alfven velocity, i.e., the amplitude of the 
magnetic fields. Therefore both linear and nonlinear calculations give almost identical results, which are AT TT (l) oc 1°, 
between 2000(~ l s ) < I < 18000. 

On small scales, I > 18000, the linear calculation shows steeper rise of the power spectra for both Cases A and B 
than the nonlinear calculation as shown in FIG. 4. There also exists exponential damping in the linear calculation on 
smaller scales. On the other hand, nonlinear calculation shows less rise and no damping in the power spectra. 

The reason why the power spectra of the nonlinear calculation have less prominent peaks on the small scales 
(large Vs) is due to the cascade decay of the magnetic fields during the turbulent regime. As shown in FIG. 2, 
the peak of the magnetic field power spectrum shifts to larger scales. The condition of this shift is determined by 
teddy = Li n t /va = 1 / H as is discussed in Sec. Ill Al Then the cutoff scale (or integral scale) has been frozen since the 
transition from the turbulent regime to the viscus regime. Accordingly the comoving wave number of the cutoff can be 
written and fei n t = aH/vA\z t ^ v as is obtained in Eq. (|13|) . On the other hand, the cutoff scale of the linear calculation 
is always determined by the diffusion condition k c — k$/vA (see Appendix). The power spectrum increases toward 
small scales as long as k < k c . Then it starts to show exponential damping at k c . In Cases A and B, fc; nt < k c at 
the recombination epoch. Therefore nonlinear calculation shows less peaks in the power spectra for these cases. The 
diffusion damping scale is determined by k c = kg/vA oc 1/B for both linear and nonlinear calculations. Therefore 
one might expect to have same damping behaviors for both linear and nonlinear calculations. However, since the 
magnetic field strength of the nonlinear calculation on the scales smaller than 1 / k mt is smaller than the one of the 
linear calculation as well as temperature anisotropics, we expect to have the smaller diffusion damping scale for the 
nonlinear calculation than the linear calculation, i.e., /£ onlinear > fe^ near . Within the range of FIG. 4, we can only 
see the diffusion (exponential) damping of the linear calculation. The nonlinear calculations for Cases A and B only 
show much milder damping due to the cascade decay during the turbulent regime below the integral scale. 

In Case C, k m t > k c at the recombination epoch. Therefore we ought to obtain very similar temperature power 
spectra for both linear and nonlinear calculations. 

VI. CONCLUSION 

In this paper we study the effect of the primordial magnetic fields on the CMB temperature and polarization 
anisotropies. In particular the nonlinear evolution of the magnetic fields and the resulting Alfven modes of the fluid 
velocities, which are the source of temperature anisotropies and polarization, are appropriately included based on the 
cosmological MHD simulation by Banerjee and Jedamzik 0, ITU . Diffusion and direct cascade processes are properly 
taken into account. 

We separate the evolution into three regimes, i.e., turbulent, diffusion, and free streaming. In the turbulent regime, 
the advective term, which is essentially nonlinear, provides the dominant contribution on the evolution of Alfven 
modes. Viscosity plays an important role in the viscous regime, while the drag term takes over the task of the viscous 
term in the free streaming regime. 
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By combining the numerical simulations of three regimes, we obtain a comprehensive evolution history of the 
magnetic fields and Alfven modes. We find the relation between the integral (or coherent) scale on which the magnetic 
energy peaks, and the maximum magnetic field strength, or equivalently, the magnetic field strength at the integral 
scale as shown in FIG. El 

We divide the evolution of the magnetic fields and the Alfven modes into three cases, i.e., Cases A, B, and C. For 
Case A, the integral scale is still in the viscous regime at the recombination epoch, while the integral scale is in the 
free streaming regime for Cases B and C. The direct cascade process takes place during the turbulent regime for all 
cases. In the free streaming regime the direct cascade process as well as the diffusion process works for Case C at the 
integral scale. The resultant velocity spectra of Alfven modes are quite different between three cases. 

From the velocity spectra, we calculate CMB temperature anisotropy and (B-mode) polarization spectra. We make 
comparisons between nonlinear and linear calculations and find differences on both large and small scales. On large 
scales / < 2000, both CMB anisotropy and polarization spectra have flat and blue spectra for the nonlinear and 
linear calculations, respectively. This difference is caused by the inclusion of the advective term for the nonlinear 
calculations. This difference gives stronger constraints for the nonlinear case in the intermediate angular scale using 
CMB observations such as WMAP and BOOMERANG. Note that the possible non-equipartition in the turbulent 
regime may make the difference between nonlinear and linear calculations small. 

Using WMAP and BOOMERANG, we set rough upper limits for the comoving magnetic field strength as B < 
28nGauss and the comoving integral scale as k- ln t > 1.5Mpc _1 if the equipartition is valid in the turbulent regime. If 
this assumption is violated, we still can set rough upper limits as B < 30nGauss and fc; nt > 1.2Mpc . 

For the polarization spectra, we also expect higher signals on large scales by nonlinear calculations than linear ones. 
The signal may exceed the B-mode polarization from the gravitational lenses if primordial magnetic fields exist. 

On small scales, nonlinear calculations show milder damping of temperature anisotropies and polarization than 
linear calculations in Cases A and B. We expect to have both temperature anisotropy and polarization spectra even 
beyond I > 10000. The peak values of temperature anisotropy and B-mode polarization spectra are approximately 
iOfiK and a few /LtK depending on the peak scale h n t (or the integral scale feint)- 

Note that we consider only the case with the spectral index of the magnetic field spectrum n > 3/2 in this paper, 
while the extension to the case with < n < 3/2 is straightforward. 

Various observation projects for the small scale temperature anisotropies and polarization are planed |32|. These 
observations may find the evidence of primordial magnetic fields or at least will set stringent constraints. Perhaps the 
best constraint will be provided by the really fine scale data, i.e., I > 10000. 

Finally, another and possibly stronger constraint on the field strength of putative primordial magnetic fields might 
be derived from the deflection of ultra high energy cosmic rays (UHECR). Dolag et al. [3j| simulated the evolution 
of magnetic fields using a MHD SPH code and solved the propagation of UHECR. Their seed magnetic fields were 
only 2 x 10~ 12 Gauss and they found that strong deflections of UHECR with arrival energies E = 4 x 10 19 within 
a distance of lOOMpc occur only when UHECR cross galaxy clusters. Therefore the arrival directions of UHECR 
mostly trace the source positions on the sky. On the other hand, if there exists the magnetic fields of lOnGauss as a 
seed, we expect to have very strong deflections for UHECR with arrival energies even above 10 20 eV. In near future, 
a huge amount of UHECR will be observed by the air shower experiments such as the Pierre Auger observatory |35| 
whose angular resolution is ~ 0.6 degree J36|. If the observed UHECR with 10 20 eV, which only can arrive from nearby 
sources (< lOOMpc) due to GZK cutoff [33, do trace the large scale structure, we will be able to set a very stringent 
upper limit for the primordial magnetic field strength. Note that the propagation and the deflection of UHECR are 
still open questions. Contrary to Dolag et al., Sigl et al. found larger deflection of UHECR in their simulation |34| . 
They claimed that the sources are strongly magnetized (~ /iGauss) and the deflection angle can be of order 20° up 
to 10 20 eV even if the extragalactic magnetic fields of the observer are negligible ( <§; 0.1/xGauss). If this is the case, 
however, we can still set an upper limit to the magnetic filed strength, if the observed UHECR above 10 20 eV do trace 
the large scale structure. In near future, a huge amount of UHECR will be observed by the air shower experiments 
such as the Pierre Auger observatory [35( whose angular resolution is ~ 0.6 degree l36[. If the observed UHECR 
with 10 20 eV, which only can arrive from nearby sources (< lOOMpc) due to GZK cutoff [37], do trace the large scale 
structure, we will be able to set a very stringent upper limit for the primordial magnetic field strength. 
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APPENDIX A: LINEAR PERTURBATION 



In this appendix, we review the linear perturbation analysis of the magnetic fields in the early universe. For a 
detailed discussion, see the references [T3, ■ 

We assume that the metric perturbations of the vector mode are expressed as Eq. Q58p. The stress energy tensors 
are divided into the fluid-part, T^ v , and the magnetic field part, Tg" . In the fluid part, we take into account of the 
viscosity between photons and baryons (electrons) 

T-p" = ( Pt + Pt)U»U" + Pt g^ - vH^H^Wafs, (Al) 



Here is the four-velocity and is written as 



U° = l/a, U l =u l /a, 



(A2) 



(A3) 



where u 1 are velocities. Now we are interested in the vector mode of the stress energy tensor perturbations, so that 
the velocity filed is divergence free and is decomposed as Eq. (IdT7|| . 
The magnetic part of the stress energy tensor is expressed as 
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(A4) 



For obtaining the vector mode of the magnetic part, we introduce a projection tensor onto the transverse plane as 

P^^S^-hkj/k 2 . (A5) 
Then the vector mode of Eq. (|A4() can be written as 

■^Bij = {Pinkj + Pjnki)k m TBmn/ k" . (A6) 

From the conservation law of the energy momentum tensors Tp and Tg , we obtain the Eulcr equation 

fcIL 



{■in - Vi) + (1 - 3c st )-(u z - V t ) + -— 

a a p t + pt 



{ui - Vi) 



a 4 (pt +Pt)' 



(A7) 



where dots represent the derivatives with respect to the conformal time and the comoving pressure gradient term of 
the magnetic field H is defined as 
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(A8) 



As discussed in Sec. IIVI we use which is defined as u, = Ui — Vi. We can rewrite Eq. (|A7|) in terms of Vi- When 
we evaluate Eq. i|A7|l at large scales where the viscosity can be neglected as is the case of the turbulent regime in the 
nonlinear calculations, we obtain as the following approximation in the radiation dominated epoch 



Vi 



A(R+l) Pl a A 



(A9) 



where R = 3/9b/4/9 7 . We refer the scales where Eq. (1A9|) is valid as "no diffusion scales". 

At small scales where the viscosity is dominant over the cosmological expansion, which we call as "viscous scales" , 
we can neglect the expansion term in Eq. (|A7(I and obtain 
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with employing the terminal velocity approximation. 

Matching Eqs. l|A9)l and (|AT0|) . we acquire the transition scale between the no diffusion scales and the viscous 
scales as 
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(All) 
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At the smaller scales than the photon mean free path, which we refer as "free streaming scales", baryons and 
photons are decoupled and we can no longer adopt the diffusion approximation, Eq. (|A7() . Instead, we must introduce 
the drag force term TVi/R in the baryon Euler equation as 

d r fell, 
Vi + -Vi + -Vi = j. A12 

When we neglect the cosmological expansion term and apply the terminal velocity approximation we obtain 



Aa 5 p 7 



(A13) 



The transition scale between the viscous scales and the free streaming scales is kf as discussed in Sec. IIII Bl because 
Eqs. (|A10|I and l|A13|) are identical to Eqs. lO and JT7J). 

For the comparison with the nonlinear results we calculate the power spectrum of the velocity fields. From the 
above results Vi is rewritten as 

Vi = kX(r])Ui, (A14) 



( 3r/4 / 9 7 a 4 (l + i?), k < k v , 
X(rj) = < 15/4p 7 a 3 /c 2 L m f p , k v < k < k { , 
{ 3L m f p /4p 7 a 5 , k { < k. 

The ensemble average of the velocity fields leads to 

(\v l \ 2 )=k 2 X(r 1 f(\U l \ 2 ). 

By following the procedure of ref. 0], (Ind 2 ) can be solved as 
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where q = |q| and /i = k • c{/(qk). Now, we assume that the magnetic power spectrum is the power law spectrum with 
the spectral index n as 



k 3 

Vb = —(\B c 



(A18) 



where the strength of the magnetic fields is B_9nGauss at the wave length k n . 

Assuming the power law spectrum for the magnetic fields, we can obtain the analytic approximation of (|n,| 2 ) 
fill fill |22| . Accordingly, the power spectrum of the fluid velocity leads to 



^( fc ) = |^(K| 2 )« fc2i? -9^W 2/ ( fc ) 2 ' 



(A19) 



where the mode coupling I(k) is different for the spectrum with the index n < 3/2 and n > 3/2. From the reference 
[l6[.in the case of n < 3/2, the mode coupling I(k) is approximated as 



I\k) 



1 



In the case of n > 3/2, 



I\k) 



k 



32-7T 4 3n V k n 



7 



2n 



32tt 4 15(2n-3) \k n 



2n-3 



(A20) 



(A21) 



where k c is the cutoff scale defined as k c = ks/vA [3, |9j . This cutoff of the power spectrum is caused by the dissipation 
due to the drag force. 
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For the comparison with the nonlinear results in this paper, we focus on the case with the blue spectrum n > 3/2. 
Using Eq. (|A21|I . we obtain the power spectrum of the fluid velocity as 



P v (k)~k 2 



3r? 



1 2 



4p 7 a 4 (l + i?) 



1(h) 2 oc fc 5 



(A22) 



V v (k) w k 2 



15 



4/9 7 a 3 /c 2 L mfp 



7(fc) 2 cx fc, k v < k < hi, 



(A23) 



V v {k) « fc 2 



3L 



mfp 



4p 7 a 5 



J(fc) 2 oc fc 5 , fc f < fc. 



(A24) 



From the velocity power spectrum, CMB temperature anisotropy and polarization spectra can be calculated. Sub- 
ramanian and Barrow 0, |22] obtained the following results with using the small angle approximation. 
In the case of fccr <C 1, 



fc=Z/(fO->)LSs) 



(A25) 



2tt 



2 V 3a 



fe = V(')0-'?LSs) 



And for the case of ka ^> 1, 



(A26) 



at tt = T 



l{l + l)C(l) 
2tt 



To 



7T 1 / 2 

2ka 



Pv(k) 



fc=(/( r /0-»)LSs) 



(A27) 



j(j + i)c(0 



AT BB = T \/ " v " ^To^/— P„(A) 



'tt 1 / 2 /fcL mfp 



2tt 



2fccr V 3a 



fc=i/(r/ -»?LSs) 



(A28) 



The temperature anisotropy spectrum induced by the magnetic fields with spectral index n > 3/2 in the low 
multipoles is given by substituting Eq. I|A22() into Eq. (|A25(1 and the one in the high multipoles is given by Eqs. 
(jA"23|) and (£53) into Eq. ljA"27)l as 



AT ' ' * :,4 /Ji M f ) /(fcO [fiK] ck ; 5 / 2 , I < l v , 



(A29) 



AT TT w 13 



•0 B 2 _ 9 (-^) 1 1(h) [fjK] ex fc°, ; v < I < k, 



(A30) 



AT TT « 0.4 £ 2 



20000 



1/2 



/(Al) [/tK] oc Z 2 , i f < I < l c , 



(A31) 



/(fc/) = I(k)\k=l/(r, -T, hSS ), 



(A32) 



where ? v = k v (r/ - ry LS s) and l c = k c (i] - ?7lss)- 

The B-mode polarization spectrum induced by the magnetic fields with spectral index n > 3/2 is acquired by 
substituting Eq. (|A22|) into Eq. I|A26|) and substituting Eq. <|A23|) into Eq. (|A28|) as 



AT BB » 0.04 B 2 



I 

1000 



1(h) \jjlK] oc f' 2 , I < l v , 



(A33) 



27 



AT BB w 0.12 5^ 



-9 




I \- 1/2 

—J /(A,) [mK] cx /, * v < I < k. 



(A34) 
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